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ABSTRACT 

The probability distributions (PDFs) of the differences of any physical vari- 
able in the intermittent, turbulent interplanetary medium are scale dependent. 
Strong non-Gaussianity of solar wind fluctuations applies for short time-lag space- 
craft observations, corresponding to small-scale spatial separations, whereas for 
large scales the differences turn into a Gaussian normal distribution. These char- 
acteristics were hitherto described in the context of the log-normal, the Castaing 
distribution or the shell model. On the other hand, a possible explanation for 
nonlocality in turbulence is offered within the context of nonextensive entropy 
generalization by a recently introduced bi-kappa distribution, generating through 
a convolution of a negative-kappa core and positive-kappa halo pronounced non- 
Gaussian structures. The PDFs of solar wind scalar field differences are computed 
from WIND and ACE data for different time lags and compared with the charac- 
teristics of the theoretical bi-kappa functional, well representing the overall scale 
dependence of the spatial solar wind intermittency. The observed PDF charac- 
teristics for increased spatial scales are manifest in the theoretical distribution 
functional by enhancing the only tuning parameter k, measuring the degree of 
nonextensivity where the large-scale Gaussian is approached for k — > 00. The 
nonextensive approach assures for experimental studies of solar wind intermit- 
tency independence from influence of a priori model assumptions. It is argued 
that the intermittency of the turbulent fluctuations should be related physically 
to the nonextensive character of the interplanetary medium counting for nonlocal 
interactions via the entropy generalization. 

Subject headings: interplanetary medium — plasmas 
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1. Introduction 

Spatial intermittency in fully developed turbulence is an established feature of astro- 
physical plasma fluctuations and is in particular manifest in the interplanetary medium by 
in situ space probe observations. Hence, the analysis of probability distribution functions 
(PDFs) or the computation of higher order moments are of considerable interest to study 
multi-scale statistical properties of the solar wind plasma and interplanetary magnetic field 
fluctuations. Classical statistical theory provides the phase space distribution from ideal 
MHD invariants, introduced by Matthaeus & Goldstein (1982) to characterize spectral prop- 
erties of incompressible solar wind fluctuations. The existence of multifractal structures in 
the fluctuations of the magnetic field strength and the plasma parameters as temperature, 
density and velocity of the interplanetary magnetofluid was confirmed (Burlaga 1991, 1992), 
followed by analyses of non-Gaussian PDFs, related intermittency and the fractal scaling of 
the solar wind MHD fluctuations by Marsch & Tu (1994, 1997). The fractal/ multifractal 
nature of certain solar wind and magnetospheric key parameters, however, remains to be a 
subject of recent investigations (Hnat et al. 2003). Interplanetary observations were used to 
establish the differences between fluid and MHD turbulence (Carbone, Veltri & Bruno 1995) 
and the Castaing distribution (Castaing, Gagne & Hopfinger 1990; Castaing & Dubrulle 
1995; Castaing 1996) was widely considered in favor of models of the non-Gaussian cascade 
character of the intermittency in turbulent flows (Consolini & De Michelis 1998; Guowei 
et al. 1998; Sorriso-Valvo, Carbone & Veltri 1999; Schmitt & Marsan 2001). In addition 
to other solar wind parameters as velocity and density the non-Gaussian and intermittent 
properties of the interplanetary magnetic field fluctuations have shown to represent a new 
set of intermittency parameters controlling the energy input rate into the magnetosphere 
(Voros et al. 2002a; Voros & Jankovicova 2002b). Multi-scale intermittency and anisotropy 
effects appear to be decisive factors in near-Earth magnetotail dynamics as well (Voros 2003, 
2004). In the solar wind, based on WIND, ACE and Voyager observations Burlaga et al. 
(2002); Burlaga (2002) (and references therein) provide a comprehensive study of multiscale 
statistical properties of speed fluctuations at 1AU and 60 AU, respectively, verifying that 
the PDFs at small scales have a shape that is characteristic for intermittent turbulence with 
relatively large skewness and kurtosis and approach a Gaussian for large scales (see also 
Sorriso-Valvo, Carbone & Veltri (1999)). 

On observational grounds the evolution of energy spectra clearly demonstrates the oc- 
currence of nonlinear interactions in the solar wind fluctuations. These interactions include 
turbulent cascades (Tu & Marsch 1995) as well as nonlocal interactions between separated 
scales (Chang 1999). 

The classical Boltzmann-Gibbs extensive thermo-statistics constitutes a powerful tool 
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when microscopic interactions and memory are short ranged and the environment is a Eu- 
clidean space-time, a continuous and differentiable manifold. However, in the present sit- 
uation we are dealing with astrophysical systems, generally subject to spatial or temporal 
nonlocal interactions evolving in a non-euclidean, e.g. multi-fractal space-time that makes 
their behavior nonextensive. A generalization of the Boltzmann-Gibbs-Shannon entropy for 
a statistical equilibrium was introduced by Tsallis (1988), preserving the usual properties 
of positivity, equiprobability and irreversibility, but suitably extends the standard extensiv- 
ity or additivity to nonextensivity. With regard to plasma astrophysics we emphasize that 
the pseudo-additive entropy concept generates power-law distributions where a parameter 
q characterizes the fraction of the non-thermal particle components. On the other hand, 
space physicists recognized four decades ago from in situ satellite observations frequently 
the existence of high energy tail velocity space distributions of electrons and ions in space 
plasmas. The family of k— distributions, a power law in particle speed, was introduced on 
phenomenological grounds to model these suprathermal structures, where k assumes only 
positive values. In continuation, significant progress was provided by Treumann (1999a,b) 
who developed a kinetic theory demonstrating that power-law velocity space distributions 
are a particular thermodynamic equilibrium state. Relating the parameters q and k by an 
elementary transformation Leubner (2002) provided the missing link between nonextensive 
distributions and k— functions favored in space plasma physics, leading to the required the- 
oretical justification for the use of k— distributions from fundamental physics. Since the 
parameter k, a measure of the degree of nonextensivity of the system, is not restricted to 
positive values in the nonextensive context, the commonly observed core-halo twin character 
of the interplanetary electron and ion velocity space distributions was verified theoretically 
upon generalization to a bi-kappa distribution, subject to a less pronounced core along with 
extended tails as compared to a Maxwellian (Leubner 2004a,b). 

Recently, the PDF of the Tsallis ensemble was linked to the analysis of fully developed 
turbulence providing a relation between nonextensivity and intermittency, which manifests 
the multifractality in the distribution of eddies (Arimitsu & Arimitsu 2000a, 2001). It was 
shown that the value of the index q of nonextensive statistics is related to the extremes of the 
multifractal spectrum (Lyra & Tsallis 1998) and that the multifractal spectrum correspond- 
ing to the Tsallis statistics is determined self-consistently (Arimitsu & Arimitsu 2000b). 
Moreover, the context of generalized thermo-statistics provides analytical formulas for PDFs 
of distance dependent velocity differences, implying the cascade like structure of the turbu- 
lent dynamics and an interpretation of the parameter k = l/(q — 1) (Beck 2000). Links 
between the nonextensive parameter q and the corresponding parameters of the log-normal, 
multifractal and random-/? models were derived (Shivamoggi & Beck 2003). 

We relate, in the following, nonlocality in turbulent flows to nonextensive systems and 
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demonstrate in the context of entropy generalization the consistency of the theoretical bi- 
kappa distribution (Leubner 2004a) with observed, scale dependent PDFs of characteristic 
key variables in the intermittent, turbulent interplanetary medium, where k appears as the 
only tuning parameter between the scales. 



2. Theory 

The nonextensive entropy generalization for the classical thermo-statistics proposed by 
Tsallis (1988) takes the form 

1 - yv 

where pi is the probability of the i th microstate, k B is Boltzmann's constant and q is a 
parameter quantifying the degree of nonextensivity of the system, commonly referred to as 
the entropic index. A crucial property of this entropy is the pseudo-additivity such that 



S q (A + B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B) (2) 

for given subsystems A and B in the sense of factorizability of the microstate probabili- 
ties. Hence, nonlocality or long-range interactions are introduced by the multiplicative term 
accounting for correlations between the subsystems. In order to link the Tsallis q— statistics 
to the family of k— distributions applied in space and astrophysical plasma modeling we 
perform the transformation 

(3) 



l-q 

to Eq. (1) yielding the generalized entropy of the form (Leubner 2002) 

S K = Kk B {Y,p]~ 1/K - 1) (4) 

where k — oo corresponds to q — 1 and represents the extensive limit of statisti- 
cal independency. Consequently the interaction term in Eq. (2) cancels and the classical 
Boltzmann-Gibbs-Shannon entropy is recovered as 



S q = -ks^Pilnpi 



(5) 
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Eq. (4) applies to systems subject to nonlocal interactions and systems evolving in a 
non-Euclidean and multifractal space-time. A further generalization of Eq. (2) for complex 
systems composed of an arbitrary number of mutually correlated systems is provided by 
Milovanov & Zelenyi (2000) where appropriate higher order terms appear in Eq. (2). Once 
the entropy is known the corresponding probability distributions are available. 

In Maxwell's derivation the velocity components of the distribution /(v) are uncor- 
rected where In/ can be expressed as a sum of the logarithms of the one dimensional 
distribution functions. In nonextensive systems one needs to introduce correlations between 
the components, which is conveniently done by replacing the logarithm functions by a power 
law. Upon transforming the "q — exp" and "q — log" functions introduced in nonextensive 
statistics (Silva, Plastino, & Lima 1998) into the k— equivalents e K (f) = [1 + / '/ 1 k] r and 
l n K,f — K-lf 1 ^ — 1] the generalized nonextensive solution for the one dimensional distribution 
in k notation reads 



— K 

(6) 

Hence, the exponential probability function of the Maxwellian gas of an uncorrelated 
ensemble of particles is replaced by the characteristics of a power law. Note also that 
the distribution (6) can be derived entirely using general methods of statistical mechanics 
(Milovanov & Zelenyi 2000). Performing the normalization to the particle density iV yields 



fh(v) = A h 



1 \v l 

1 + —2 



N 1 ^ 

h ~ 7^7^r[ K -i/2] (7) 

and we assign to fh the notation 'halo' distribution. v t = ^/iksT jm is the thermal 
speed of a Maxwellian, where T and m are the temperature and mass, respectively, of the 
species considered. The case k = oo recovers the Maxwellian equilibrium distribution, see 
Fig. 1 and the discussion below. 

We note that contrary to Eq. (6), derived from the nonextensive entropy generalization 
Eq. (1), the functional form of the k- distributions conventionally used in plasma physics 
is powered by — (k + 1) instead of — k (Leubner 2002, 2004a). Consequently, the nonex- 
tensive distribution (6) provides more pronounced tails as compared to the conventional 
k — distribution for the same values of k. Furthermore, in all space plasma applications the 
spectral index k was hitherto assumed to be restricted to positive values, requiring in the 
derived one-dimensional nonextensive representation | < k < oo equivalent to — 1 < q < 1 
(for q < — 1 corresponding to k < 1/2 the distribution (6) cannot be normalized). 
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Since the nonextensive formalism is not restricted to — 1 < q < 1 we generalize to 
values q > 1 providing the full range — oo < k < oo for the spectral index k. Hence, upon 
incorporating the sign into the equation the one dimensional distribution for negative definite 
k can be written in the form 



fc(v) 



lv 

1 2 

K Vf 



21 K 



(8) 



where the normalization constant reads 



A = N 1 F ^ + 3/2] 

C V*V t y/t r[K+l] W 

We assign the notation 'core' distribution to f c , a function subject to a thermal cutoff 
at the maximum allowed velocity v max = \/Tzv t , which limits also any required integrations. 
Since the sign is implemented into Eq. (8), k has to be understood as a positive quantity 
wherever it appears. 

The nonextensive distribution-functionals (6) and (8), along with the proper normal- 
ization from the entropy generalization in Eq. (1), provide access to nonlocal effects, where 
k measures the degree of nonextensivity within the system. Fig. 1 demonstrates schemati- 
cally the non-thermal behavior of the suprathermal halo component fh and the reduced core 
distribution f c , subject to finite support in velocity space and exhibiting a thermal cutoff at 
v = Vt\J~K where both functions approach the same Maxwellian as k — > oo. 

With regard to a physical elucidation we refer to recent developments indicating that the 
parameter k (or equivalently q via the transformation 1/(1 — q) — k, as outlined above), as 
degree of nonextensivity in the Tsallis context, finds at least three more interpretations. In 
particular, Chavanis (2003) illuminates the Tsallis q-entropy approach in view of incomplete 
mixing in complex systems, clarifying that it is related to a particular Boltzmann H-function, 
where the parameter q measures the degree of mixing. Furthermore, the PDFs of the Tsallis 
ensemble are linked to the analysis of fully developed turbulence (Arimitsu & Arimitsu 
2000a, 2001), where the nonextensive parameter q is in particular related to nonlocality 
in the system and measures the scale dependence of the PDFs of the dynamical physical 
variables considered. Finally, the parameter q may find also a possible alternative physical 
interpretation in terms of the heat capacity of a medium. A system with n < was suggested 
to represent a heat bath with finite positive heat capacity and vice versa, for k > the heat 
bath has negative finite heat capacity. For a discussion of this property with regard to entropy 
generalization see Almeida (2001). This outline indicates that the conceptual background 
for the description of nonextensive systems requires the development of unifying concepts. 
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It should be emphasized that, besides a brief comment (Leubner 2002), the extension to 
negative k— values as formulated in Eq. (8) was just recently recognized as being physically 
relevant (Leubner 2004a,b) since research focused hitherto on the suprathermal particle 
populations, in particular, well represented by a halo distribution. Let us point out that 
the thermal cutoff of the functional Eq. (8) provides apparently a clear definition of a core 
population, whereas a 'core' of the function fh is not definable and just in some arbitrary way 
related to the central part of a Maxwellian. Besides a clearly defined core-halo transition any 
unique and physically relevant nonextensive PDF must be subject to three further conditions. 
We require that (a) the distribution approaches one and the same Maxwellian as k — > oo, 
(b) a unique, global distribution must be definable by one single density and temperature 
and (c) upon variation of the coupling parameter k particle conservation and adiabatic 
evolution should be considered, such that a redistribution in a box (a source free environment) 
can be performed. Subject to these constraints the appropriate mathematical functional, 
representing observed core-halo structures in nonextensive astrophysical environments, is 
available from the elementary combination F c h = B c h(fh + fc)- 111 this context the full 
velocity space distribution compatible with nonextensive entropy generalization and obeying 
the above constraints reads 



to which we assign the notation 'bi-kappa distribution' and where B c h = N/ (ir 1 ^ 2 Vt)G(K,) 
and the function G(k) is defined by 



from the normalization. The function G(k) is subject to a particular weak k— dependence 
where G(k) ~ 1/2, see Leubner (2004a) for a graphical illustration and discussion. Hence, 
the resulting factor 1/2 in the normalization constant B c h = N/(2tc 1 ^ 2 v t ) reflects consis- 
tently the superposition of the two counter-organizing contributions in Eq. (10). Upon 
re-distribution of particles as a consequence of nonlocal interactions, where k act as measure 
of nonlocality in the system, a fraction of core-particles is removed and available to generate 
a pronounced non-Maxwellian tail population counting for zero net particle budget. Hence, 
the specific nonextensive velocity space distribution (10) provides a unique mathematical 
core-halo representation subject to the non-trivial property of particle conservation in the 
sense that the normalization is independent of the parameter k (Leubner 2004a,b). This 
specific feature allows also to transform the one-dimensional bi-kappa velocity distribution 




(10) 




(11) 
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of Maxwell's particle context into the mathematical form of a Gaussian of constant vari- 
ance, a PDF applicable to any characteristic physical variable suitable to analyze statistical 
properties of astrophysical systems, see below. Contrary to the conventional k— distribution 
models constraint by a positive k— value, the functional form (10) permits in a closed system 
a re-distribution of energy from a Maxwellian into highly non-thermal core-halo structures 
where no support from the environment is required. This particular physical feature is only 
available if k assumes the same value for both, core and halo fraction and is manifest simul- 
taneously in a core-width reduction and distinct tail formation as k decreases. For k — oo 
the function G(k) = 1/2 and the power laws in the brackets of the right hand side of Eq. 
(10) turn each into the same Maxwellian exponential. 

Finally, besides particle conservation also adiabatic behavior is known as typical char- 
acteristic from a variety of astrophysical plasmas studies (Baumjohann & Paschmann 1989; 
Totten, Freemann, & Arya 1995). The effect of long-range interactions on the entropy depen- 
dence S(k) in a system obeying a bi-kappa PDF (10) can be studied after a decomposition 
into the core and halo fractional entropy contributions (Leubner 2004a). Nonlocality in the 
system generates, as consequence of the pronounced tails, an entropy increase with regard 
to the halo PDF that is accompanied by a counteracting partner, the core entropy con- 
tribution originating from the cutoff PDF f c , a higher organized state of reduced entropy 
as compared to the Maxwellian equilibrium entropy. In particular, the entropy of the full 
bi-kappa distribution Eq. (10) follows the Boltzmann-Gibbs-Shannon extensive behavior 
down to k ~ 10 turning after a transition into the nonextensive profile of fh alone. Since 
dS > as k decreases from k — oo to low values any generation of highly non-Maxwellian 
PDF structures due to nonlocality in multifractal systems obeys consistently the second law. 

Intermittency in astrophysical plasma turbulence can be analyzed by comparing ob- 
served PDFs of fluctuations of the velocity field magnitude in the solar wind with the global 
theoretical core-halo representation (10), accounting for nonlocality through entropy gener- 
alization. On the other hand, on small scales also the PDFs of interplanetary density and 
magnetic field fluctuations are observed to obey the same leptokurtic core non-Gaussianity 
along with pronounced tails or halos, requiring a generalization of Eq. (10). 

We introduce a unique, global PDF provided in the context of nonextensive thermo- 
statistics for any physical variable, where nonlocality is associated to pseudo-additive entropy 
generalization. Upon normalizing the one-dimensional bi-kappa particle distribution (10) to 
unity and assigning a unique distribution variance a to the thermal spread v t , the 'Maxwellian 
form' of the bi-kappa distribution (10) transforms with G(k) = 1/2 to a 'Gaussian form' of 
a global bi-kappa PDF as 
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P cH (SX;k) 
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(12) 



Subject to a single constant variance a this one-parameter PDF is applicable to the 
differences of the fluctuations 5X(t) = X(t + r) — X(t) of any physical variable X in the 
astrophysical system considered, where X(t) denotes any characteristic solar wind variable 
at time t and r is the time lag. Here k assumes a clear physical interpretation defining the 
degree of nonextensivity or nonlocality in the system, thus being a measure of long-range 
interactions. As k — ► oo the bi-kappa distribution P ch (SX; n) approaches a single Gaussian. 
In the following we demonstrate the relevance of the nonextensive context by analyzing the 
observed scale dependence of the PDFs of the differences of velocity, density and magnetic 
field variables in the intermittent, turbulent interplanetary medium based on the theoretical 
global bi-kappa PDF (12). 



3. Data Evaluation 

For a proper estimation of empirical PDFs of interplanetary plasma and magnetic field 
data long sets are needed. Due to non-stationarity in driving dissipation conditions, however, 
the data sets have to be chosen carefully. For the present approach solar wind velocity data 
are analyzed with a time resolution of 92 s and the interplanetary magnetic field strength 
with a time resolution of 60 s, both available from the WIND spacecraft. WIND has a 
complicated trajectory that also crosses the magnetosphere from time to time. Our analysis 
is based on the entire year 1997 data set when WIND was in the GSE positions: 150-230 Re- 
This period is compared with a shorter time period from November 18 to December 10, 1998 
when magnetic field measurements with time resolution of 16 s were available from the ACE 
spacecraft. The ACE satellite is continuously monitoring the solar wind at the LI point. 
Data gaps were interpolated in all cases and the resulting PDFs are fitted by the theoretical 
nonextensive distributions. Moreover, we provide a comparison of the WIND magnetic field 
analysis with the corresponding PDFs obtained from ACE magnetic field measurements. In 
particular, from the 1997 WIND velocity and magnetic field data 5 day periods of low speed 
subintervals were selected such that during each subinterval the bulk velocity was less than 
450 km/s. 

No velocity observations were available from ACE during the interval of January 10- 
29, 1998. However, at this time WIND was upstream, close to the ACE position. WIND 
measurements provide clear evidence that between November 18 and December 10, 1998, the 
bulk speed was less than 450 km/s in about 92 % of this period. This interval was treated in 
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one block as a whole since also the magnetic field intermittency was low, where any increase 
of the latter can be associated with high speed streams and increased geomagnetic response 
(Voros et al. 2002a). 

For each data set the plasma and magentic field increments were calculated at a given 
time lag r, equivalent to a particular spatial scale by 5X(t) = X(t + r) — X(t) where each 
variable within each subinterval is normalized to the standard deviation. For each 5 day long 
realization of the low speed flow the empirical probability distribution function (histogram) 
was then computed. SX(r) is binned into iV equal spaced boxes and the number of elements 
in each box was computed where the robustness of the histograms against N is tested. 5X(t) 
represents characteristic fluctuations at the time scale r or equivalently across eddies of size 
I ~ V * r. In the case of single spacecraft measurements a similarity in statistical properties 
of temporal or spatial fluctuations is supposed. This is true in particular for the solar wind 
(Marsch & Tu 1997) where the spatial fluctuations on a scale I pass over the spacecraft 
faster than they typically fluctuate in time. Hence, by changing r it is possible to analyze 
the statistical features of fluctuations in different time scales, which roughly correspond to 
those statistical characteristics across turbulent eddies of size I ~ V * r. 

At large scales, e.g. r = 24 hrs. and I ~ 5000i?£ (supposing V = 350 km/s) the 
PDF's are close to the Gaussian distribution, because the correlation between the separated 
points vanishes, whereas at small scales, e.g. r = 0.2 hrs. and / ~ AORe, the PDF's clearly 
have non-Gaussian shapes, see section 4. These values were also estimated by Sorriso-Valvo, 
Carbone & Veltri (1999). 

4. Results and Discussion 

Core-halo structures are a persistent feature of a variety of different astrophysical sys- 
tems. In particular, a manifestation of those formations is found in the thermo-statistical 
properties of interplanetary plasmas (e.g. Marsch et al. (1982); Pierrard, Maksimovic, & 
Lemaire (1999); Maksimovic, Gary, & Skoug (2000)), in gravitational equilibrium states of 
stellar system (Chavanis & Sommeria 1998) or in the distribution of dark matter in galaxies 
and beyond, see e.g. Burkert (1995); Salucci & Burkert (2000). We focus in this study 
on astrophysical plasma turbulence, specifically on PDFs of fluctuations in the intermittent 
solar wind turbulence, safely available from in situ satellite observations. 

The Wind and ACE solar wind data analysis unambiguously manifests that the PDFs 
of large scale velocity, density and magnetic field fluctuation differences are well represented 
by a Gaussian, turning simultaneously into leptokurtic peaked distributions of strong non- 
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Gaussianity in the center along with a pronounced tail structure on small scales. In par- 
ticular, the PDFs of large-scale magnetic field fluctuations, not related to the increment 
field are known to be subject to relatively small deviations from the Gaussian statistics and 
are well fitted by the Castaing distribution, a convolution of Gaussians with variances dis- 
tributed according to a log-normal distribution (Castaing, Gagne & Hopfinger 1990; Padhye 
et al. 2001). Independent of the physical situation considered, the Castaing distribution 
provides a multi-parameter description of observed PDFs, plausible in this case, since the 
large scale fluctuations of the interplanetary magnetic field are generated by a variety of 
discrete coronal sources. If individual coronal sources evoke Gaussian distributed magnetic 
fields, the net magnetic fluctuations can be modeled by their superposition with a spread of 
the corresponding variances. Here the Castaing model is not related to any energy cascading 
processes. 

Contrary, small-scale fluctuations of the increment field are associated with local inter- 
mittent flows where fluctuations are concentrated in limited space volumes. Consequently, 
the PDFs are scale dependent and intermittency generates long-tailed distributions. It is 
customary to use n-th order absolute powers of the plasma variables and magnetic field 
increments (n-th order structure functions, Pagel & Balogh (2001); Marsch & Tu (1997)) 
allowing to investigate the multi-scale scaling features of fluctuations subject to long-tailed 
PDFs. Direct studies of observed PDFs of the increment fields 5X(t) = X(t + t) — X(t) 
for any characteristic solar wind variable at time t and time lag r revealed departures from 
a Gaussian distribution over multiple scales (Sorriso-Valvo, Carbone & Veltri 1999) and an 
increase of intermittency towards small scales (Marsch & Tu 1994). The PDFs are also 
found to be leptokurtic which indicates the turbulent character of the underlying fluctua- 
tions. Sorriso-Valvo, Carbone & Veltri (1999) have shown that the non-Gaussian behavior 
of small-scale velocity and magnetic field fluctuations in the solar wind can also be described 
well by a Castaing distribution, where the individual sources of Gaussian fluctuations appear 
as a fragmentation process through scales in turbulent cascades. 

A more general fragmentation process based on the extreme deviation theory was pro- 
posed by Frisch & Sornette (1997). In this model the increment fields are due to the product 
of fragments derived from structures at larger scales and the PDF is a stretched exponential 
function P se (5X) ~ exp[— /3|5X| M ], where (5 and \i are the tunable parameters, governing the 
fragmentation process. Sorriso-Valvo et al. (2001) have shown that both, the Castaing and 
the stretched exponential model can describe intermittency in the solar wind equally well. 
Fig. 2 provides a comparison of the stretched exponential model, the widely used log-normal 
model (Frisch 1996) Pi n (5X) ~ 1 / (a5X)exp[—ln(5X) 2 /a 2 ], where a is the variance of the 
distribution and m is the scale parameter, both compared to a bi-kappa fit to ACE mag- 
netic field amplitude data between November 18 and December 10, 1998. The dimensionless 
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t = 14 is multiplied by the time resolution of 16 sec. to generate an effective time-lag and 
the log-normal model restricts to positive values of the differences. All three models pro- 
vide excellent fits to the data indicated by the error bars for large differences, whereas the 
stretched exponential model slightly deviates from experimental curve at small SB. 

In Fig. 3 we focus on the bi-kappa model showing for the same date set that the scale 
dependence of the PDF in the solar wind can be represented accurately via the only tuning 
parameter k, which measures the degree of nonextensivity or coupling within the astrophysi- 
cal plasma environment. The nonextensive context, generates a precise representation for the 
observed PDFs characterizing the intermittency of the small scale fluctuations. For compar- 
ison also a Gaussian and the conventional k— function (Leubner 2002), subject to the same 
k — value and not able to reproduce the structure of the PDF of small scale fluctuations, are 
provided. Here and in the following figures the dotted lines indicate the standard deviation. 

Based on WIND velocity field magnitude data Fig. 4 presents an analysis of the scale 
dependence of interplanetary PDFs of the velocity field magnitude indicating in three plots 
from left to right the decreasing intermittency with increasing spatial separation scale or 
time-lags, where now the observational uncertainty is indicated by the standard deviation 
(dotted lines). Constraint by a time resolution of 92 sec. the two point time separation is 
t x 92sec where r assumes the values 10, 70 and 900 from left to right. The solid lines repre- 
sent best fits to the observed PDFs with nonextensive distributions where the corresponding 
k is determined as k — 2, 3.5 and oo. The strong non-Gaussian character of the leptokurtic 
PDFs (left panel) exhibiting pronounced tails associated with solar wind turbulence and 
intermittency in small scale fluctuations finds again an accurate analytical fit and hence a 
physical background in the nonextensive representation, accounting for nonlocality in tur- 
bulence. The non-Gaussian structure is somewhat softened for enhanced r = 70 (central 
panel) but again precisely modeled within the pseudo-additive entropy context turning into 
the Gaussian shape of large scale fluctuations, independent of the increment field. 

Fig. 5 and Fig. 6 provide the corresponding nonextensive analysis of the scale de- 
pendence of the density and magnetic field fluctuations obtained from WIND data. Ev- 
idently, the scale dependent characteristics of the observed PDFs of the increment fields 
SX(t) = X(t + r) — X(t) for all solar wind variables evolve simultaneously on small scales 
approaching independency of the increment field in the large scale Gaussian. With separa- 
tion scales of r = 10, 70, 900 the corresponding evolution of the PDFs of observed density and 
magnetic field fluctuations are best represented by the same values of k — 2, 3.5, oo counting 
for nonlocal interactions in the nonextensive theoretical approach as for the velocity field 
magnitude. 

Finally, Fig. 7 demonstrates from ACE magnetic field observations that the results of the 
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structure analysis regarding the scale dependence of the increment field PDFs of physical 
variables in the interplanetary medium performed by WIND data is not restricted to the 
observations of this particular spacecraft. Highly accurately, the overall scale dependence 
appears as a general characteristic of quiet astrophysical plasma environments. In particular, 
the scale dependent PDF evolution of magnetic field fluctuations is subject to a two point 
separation scale of r = 100, 2000, 10000 with a time resolution of 16 sec. Best fits of the 
bi-kappa distribution are obtained for k = 1.8, 3, oo, measuring the degree of nonextensivity, 
or coupling through nonlocal interactions, respectively. 

Figs. 4, 5 and 6 demonstrate clearly via the PDFs of the density and the velocity and 
magnetic field magnitudes that the scale dependence of the fluctuations of different charac- 
teristic interplanetary variables exhibits the same behavior or is at least similar. The analysis 
is based on slow speed solar wind PDFs indicating a universal scaling dependency between 
N,v and B intermittency within the experimental uncertainties. This strong correlation im- 
plies that the scale dependencies of all physical variables are coupled where the solar wind 
Alfvenic fluctuations provide a physical basis of the velocity and magnetic field correlations. 
Furthermore, the proportionality between density fluctuations and the magnetic field and 
velocity fluctuations is already maintained in the solar wind by the presence of weak spatial 
gradients (Spangler 2003). 

k, — distributions reproduce the Maxwell-Boltzmann distribution for k — > oo, a situation 
identifying Kj dbS dull ordering parameter that acounts for correlations within the system. 
Quasistationary, highly correlated turbulent conditions characterized by kappa distributions 
represent stationary states far from equilibrium where a generalization of the Boltzmann- 
Shannon entropy, as measure of the level of organization or intermittency, applies (Goldstein 
& Lebowitz 2004; Treumann 2004). Physically this can be understood considering a system 
at a certain nonlinear stage where turbulence may reach a state of high energy level that is 
balanced by turbulent dissipation. In this environment equilibrium statistics can be extended 
to dissipative systems, approaching a stationary state beyond thermal equilibrium (Gotoh 
& Kraichnan 2002). However, no final theory exists for the problem of nonequilibrium 
statistics. Since turbulence is driven in the solar wind by velocity shears we choose for the 
data analysis intervals of low speed solar wind with limits in velocity space, where the driving 
and dissipation conditions do not change significantly, maintaining therefore the dynamical 
equilibrium condition. 

MHD turbulence in the solar wind is generally in a state in which propagating waves and 
non-propagating fluctuations coexist (Shebalin 1983). Moreover, typical processes include 
fluctuation induced multiscale interactions of coherent structures for which a local interaction 
in Fourier space is not required and interactions of opposite traveling waves are followed by 
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energy cascades in Fourier space. Physically, at least a part of the nonlocal forces and 
turbulent intermittency results from a source exhibiting nonlocal interactions in the Fourier 
space. We proposed a model for describing those nonlocal interactions on the basis of 
entropy generalization. However, alternative models accounting for nonlocality, responsible 
for intermittency in general, are available e.g. (Tsinober 1990; Chang 1999; Laval & Dubrulle 
2001; Falkovich et al. 2001) demonstrating that the origin of intermittency due to nonlocal 
interactions, in particular in astrophysical plasmas as the solar wind, remains open. 

Summarizing, a bi-kappa distribution family turns out theoretically as consequence of 
the entropy generalization in nonextensive thermo-statistics and was successfully applied in 
velocity space to interplanetary core-halo electron distributions as well as double humped 
proton structures. Upon transforming to the mathematical 'Gaussian form' the particu- 
lar global, on-parameter bi-kappa PDF of constant variance provides also access to the 
pronounced core-halo distributions observed in astrophysical plasma turbulence. A redistri- 
bution of a Gaussian on large scales can be performed physically in a closed environment 
turning into highly non-Gaussian leptokurtic and long-tailed structures. Pseudo-additive 
entropy generalization provides the required physical interpretation of the only tuning pa- 
rameter k in terms of the degree of nonextensivity of the system as a measure of nonlocality 
or coupling. The scale dependence of the PDFs of the differences of characteristic physical 
variables in the intermittent turbulent interplanetary medium as density, velocity and mag- 
netic field magnitude is accurately represented by the nonextensive bi-kappa functional. It 
is argued that the leptokurtic, long-tailed non-Gaussian core-halo PDFs characterizing the 
intermittency of the turbulent fluctuations can be related to the nonextensive character of 
the interplanetary medium via the entropy generalization. 
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Fig. 1. — A schematic decomposition of the characteristics of the nonextensive bi- kappa 
distribution family Eq. (10): with k = 3 the outermost curve (left panel) and innermost 
curve (right panel) correspond to the halo fh and core f c distribution fraction. For increasing 
k, — values both sets of curves merge at the same Maxwellian limit indicated as bold line, fh 
from outside and f c from inside. 
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Fig. 2. — A comparison of a stretched exponential (circles, (3 = 1.95, fi = 0.71) and a log- 
normal (triangles, a = 0.9, m = 0.75) distribution with the nonextensive bi-kappa function 
(solid line, k = 1.8) and ACE magnetic field fluctuation differences (error bars). 
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Fig. 3. — A nonextensive bi-kappa fit (solid line) with k = 1.8 of an observed PDF (dashed 
line) obtained from ACE magnetic field amplitude data. A Gaussian (thick dashed-dotted 
line) and a conventional k— distribution (thin dashed-dotted line) are provided for compari- 
son. The dotted lines indicate the standard deviation. 
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Fig. 4. — Left panel: The PDF of the increments of observed WIND velocity field magnitude 
fluctuations for r = 10 and a resolution of 92 sec. as compared to the bi-kappa function with 
k — 2. Based on the same data the central panel provides the characteristics for increase 
t = 70 where k assumes a value of 3.5 for the best representation. The PDF of large scale 
velocity fluctuations are well modeled by a Gaussian with k = oo, right panel. 
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Fig. 5. — Left panel: The PDF of the increments of observed WIND density fluctuations 
for r = 10 and a resolution of 92 sec. as compared to the bi-kappa function with k — 2. 
Based on the same data the central panel provides the characteristics for increase r = 70 
where k assumes a value of 3.5 for the best representation. The PDF of large scale density 
fluctuations with r = 900 are well modeled by a Gaussian with k = oo, right panel. 



-23- 




-5 5 -5 5 -5 5 



8 B [normalized] 8 B [normalized] 8 B [normalized] 

Fig. 6. — Left panel: The PDF of the increments of observed WIND magnetic field fluctua- 
tions for r = 10 and a resolution of 92 sec. as compared to the bi-kappa function with k — 2. 
Based on the same data the central panel provides the characteristics for increase r = 70 
where k assumes a value of 3.5 for the best representation. The PDF of large scale magnetic 
field fluctuations with r = 900 are well modeled by a Gaussian with k — oo, right panel. 
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Fig. 7. — Left Panel: the PDF of the increments of observed ACE magnetic field fluctuations 
for r = 100 and a resolution of 16 sec. as compared to the bi-kappa function with k = 1.8. 
Based on the same data the central panel provides the characteristics for increased r = 2000 
where k assumes a value of 3.0 for the best representation. The PDF of large scale magnetic 
field fluctuations, r = 10000, are well modeled by a Gaussian with k — oo, right panel. 



